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ABSTRACT: The performance of the LHCb Muon system and its stability across the full 2010 data 
taking with LHC running at = 7 TeV energy is studied. The optimization of the detector setting 
and the time calibration performed with the first collisions delivered by LHC is described. Particle 
rates, measured for the wide range of luminosities and beam operation conditions experienced dur- 
ing the run, are compared with the values expected from simulation. The space and time alignment 
of the detectors, chamber efficiency, time resolution and cluster size are evaluated. The detector 
performance is found to be as expected from specifications or better. Notably the overall efficiency 
is well above the design requirements. 
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1 Introduction 

LHCb is an experiment dedicated to heavy flavour physics at the LHC pp collider. Its primary goal 
is to look for indirect evidence of new physics in CP violation and rare decays of beauty and charm 
hadrons. 

The LHCb apparatus [1] is a single-arm forward spectrometer, consisting of a series of sub- 
detectors aligned along the beam axis. A silicon-strip Vertex Locator (VELO) centered on the 
interaction point allows for precise vertex reconstruction. A dipole magnet provides the bending 
for momentum measurement. Four multi-layer stations, one placed upstream and three down- 
stream the magnet, ensure the tracking. Silicon strips are used in the upstream station (TT) and 
in the downstream inner tracker (IT) while straw tubes are used in the downstream outer tracker 
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(OT). Particle identification is provided by two ring imaging Cherenkov (RICH) detectors, by an 
electromagnetic and hadron calorimeter system (ECAL and HCAL) and by the Muon Detector. 

Muon identification is assured by a system comprising 5 detector stations, 1380 chambers of 20 
different types for a total of 122k channels [1]. The detector and its associated readout electronics 
were designed to trigger with high efficiency on high transverse momentum muons within a time 
window smaller than 25 ns to unambiguously identify the LHC bunch crossing. 

In 2009 a first setting of the detector was performed using cosmic rays [2]. In 2010 the large 
data samples made available by the first LHC pp collisions at y/s = 7 TeV, allowed for a rapid 
improvement of the detector working conditions. 

In this paper we describe the actions taken to optimize the detector performance, in particular 
the time calibration. We then report a detailed study of the detector behaviour through the wide 
range of beam operation conditions and the five orders of magnitude of luminosity spanned during 
the 2010 data taking. The detector performance in terms of space and time alignment of the de- 
tectors, chamber efficiency, time resolution, noise level and cluster size is evaluated and compared 
with simulation and design specifications. 

Since the end of the 2010 run period considered here, LHCb has taken almost two orders of 
magnitude more data, with even higher instantaneous luminosities. The detector performance has 
remained very similar. 

2 The LHCb muon system 

The muon detector system is designed to send binary information to the data acquisition (DAQ) 
and to the hardware processors of the first level muon trigger (L0MU). 

The detector is composed of five stations (M1-M5) of rectangular shape, placed along the 
beam axis as shown in Fig. 1. Each station is equipped with 276 detectors. The detectors are 
multi-wire proportional chambers (MWPCs) with the exception of the inner part of the first station 




Figure 1: (a) Side view of the LHCb Muon Detector, (b) Station layout. The four regions R1-R4 
and the single chambers are shown. 
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equipped with triple-GEM detectors [3]. Each station consists of two mechanically independent 
halves, called A and C sides, that can be horizontally moved to access the beam pipe and the 
detector chambers for installation and maintenance. Stations M2 to M5 are placed downstream 
the calorimeters and are interleaved with 80 cm thick iron absorbers. Their information is used to 
identify and trace penetrating muons both in the online and offline analysis. Station Ml is instead 
located in front of the calorimeters and it is only used in the first level hardware trigger. 

The LOMU trigger [4] performs a stand-alone muon track reconstruction requiring hits in all 
5 stations and calculates the transverse momentum pj of the tracks. Muon candidates are accepted 
if their pj is above a given threshold. The alignment is first verified in the four stations M2-M5 
searching for hits inside suitable fields of interest (FOI) projective to the interaction point. If this 
alignment is found, the hits of M2 and M3 stations are used to predict the track hit position in Ml. 
If the Ml hit nearest to the prediction is found inside a suitable FOI, this hit and the one in M2 are 
used to define the track after the magnet deflection. The direction of such a track, its impact point 
at the magnet centre and the average pp interaction point, determine the magnet kink and the pr- 
The information of Ml station, placed in front of the calorimeter material, improves the resolution 
of the pj measurement with respect to the one that could be obtained using only the 4 downstream 
stations. The Ml information is however not used by the higher level software trigger (HLT) which 
performs a matching of the muon track detected in M2-M4 with tracks reconstructed making use 
of the full spectrometer. 

The geometry of the five stations is projective. The transverse dimensions of the stations scale 
with their distance from the interaction point. The chambers are positioned to form, across the 
stations, adjacent projective towers pointing to the beam crossing position. 

The chambers are partitioned into physical channels whose size is constrained by construc- 
tional reasons, or by requirements on their electrical capacitance and rate capability that influence 
the noise level and dead time of the front end (FE) electronics. Appropriate combinations of phys- 
ical channels are performed to build up rectangular logical pads having the x and y sizes required 
to obtain the desired performance of muon trigger and offline muon identification. 

Each station is divided into four regions with increasing distance from the beam axis. The 
linear dimensions of the regions Rl, R2, R3, R4, and the size of their logical pads, scale in the 
ratio 1:2:4:8. Since the dipole magnet provides bending in the horizontal plane, the logical pad 
segmentation of muon chambers is finer in the horizontal direction x than in the vertical direction 
y, to allow a good estimate of the momentum. Stations Ml, M2 and M3, used by the trigger to 
determine the track direction and the pj of the candidate muon, have a higher x granularity than 
stations M4 and M5, whose main purpose is the identification of penetrating particles. In the inner 
region of the first station Ml, the logical pad size is 1 cm in x and 2.5 cm in y. In the other 
stations the vertical size y just scales projectively with their distance from the interaction point; the 
x granularity instead is two times finer in stations M2, M3 and two times larger in M4, M5. The 
total number of logical pads is 55296. 

The LOMU trigger requires a five-fold coincidence among all the stations, therefore the effi- 
ciency of each station must be > 99% to obtain a trigger efficiency of at least 95%, within a time 
window smaller than 25 ns in order to unambiguously identify the bunch crossing (BX). To com- 
ply with this stringent requirement, excellent time resolution and redundancy of the detector are 
needed. The desired performance is obtained with an optimized charge-collection geometry and 
using a fast gas mixture 1 . Moreover the chambers are multi-gap detectors. In stations M2 to M5 
the MWPCs consist of two coupled bi-gap detectors with two independent readouts. In station Ml, 

lr rhe gas mixture is Ar/C0 2 /CF 4 ~ 40/55/5 for MWPCs and ~ 45/15/40 for GEM chambers. 
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R2 to R4 the MWPC's have only two gas gaps with independent readout to minimize the material 
in front of the electromagnetic calorimeter. In region M1R1 two superimposed triple-GEM cham- 
bers are used. In all cases, in standard running conditions, the two independent readout layers are 
connected to a logical "or" in the FE. 

Since constructional constraints, as well as requirements on spatial resolution and rate capabil- 
ity, strongly vary in different stations and regions of the detector, different readout techniques were 
employed. In the high granularity regions M2/3 Rl/2 a double readout was adopted for the cham- 
bers: the physical channels are narrow vertical anode strips defining the x resolution and larger 
cathode pads defining the y resolution. Both signals are readout via FE channels and directly sent 
to the trigger and DAQ through optical links. Logical pads are then obtained as a logical "and" be- 
tween anode and cathode pads. In all the other stations and regions a single readout was adopted: 
the chambers are segmented into anode or into cathode pads generally smaller than the required 
logical pad. These pads are subsequently combined to build larger logical channels that are sent 
to the trigger and DAQ. In the large low resolution external regions R4, anode pads are formed 
by soldering an appropriate number of adjacent wires. In the other regions, both in MWPCs and 
in GEM chambers, cathode pads are obtained with a segmented printed circuit board. A total of 
122112 physical channels enter the front end electronics. 

The FE boards house two custom made 8-channels ASICs called CARIOCA [5] that can pro- 
cess both the negative and positive polarity signals from wires and cathode pads. Each channel 
consists of a fast low-impedance charge-sensitive preamplifier, a main amplifier- shaper with base- 
line restorer and single-threshold discriminator. The response has a peaking time of 10 ns and a 
width not exceeding 60 ns to minimize dead time. The CARIOCA output channels are routed to 
flexible logical units performing the logical "or" of a variable number of channels. Up to four 
adjacent physical pads are connected to build a logical pad. In the Ml station, where the channel 
occupancy is high, the signals from the logical pads are directly sent to the trigger and DAQ. In 
most of the other low occupancy regions, M2/3 R3/4 and M4/5 R2/3/4, several contiguous logical 
pads are further connected to build larger logical channels in the form of vertical and horizontal 
strips, with the aim of reducing the number of optical links to the trigger and DAQ. The logical 
"or" is fully performed by logic units (DIALOG [6]) sitting on the FE boards only in part of the 
detector, while is completed on special intermediate boards in regions where the logical channel 
spans more than one FE board. 

A total of 25920 logical channels are built and finally routed to the Off Detector Electronics 
(ODE) [7], where signals are tagged with the identification number of the bunch crossing (BXID) 
and sent to the trigger processors via optical links without zero suppression. In the ODE boards 
the fine time information inside the 25 ns gate, measured by a 4-bit TDC ASIC [8], is added to the 
data to be sent to the DAQ. The trasmission from the ODE to the acquisition system goes through 
the TELL1 boards [9] where the data are zero suppressed and suitably packed. 

A distributed system, based on a CAN-bus network and 600 microcontrollers accessing sin- 
gle physical channels, performs the setting of the operating conditions and the monitoring of the 
detector [10] 

3 Detector setting 

The operating conditions of the chambers were optimized through a procedure described else- 
where [11, 12] which aims at reaching a low and stable noise rate while minimizing ageing effects 
and cluster size. These conditions were satisfied by setting the thresholds to 6 noise r.m.s. and by 
choosing the minimal high voltage value allowing to obtain the required 99% efficiency. 
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Although all the MWPCs have the same gas gap of 5 mm, different HV (from 2.53 to 2.65 kV) 
were required in different regions, depending on the pad size and the readout technique. The cham- 
ber setting was kept constant during the whole 2010 run, except for region M5R2 were voltage was 
reduced by 90 V when a large increase in luminosity produced a jump of the HV trip rate. Thresh- 
olds were lowered correspondingly to maintain high efficiency at the price of a slight increase of 
the noise level. 

The triple GEM detectors of M1R1 were operated with voltages of 435/425/415 V for the first 
part of the run. Voltages were later reduced twice by 5 V per gap for safer operations following the 
increasing luminosity. 

4 Detector Operation in 2010 

During 2010 the detector was operated to acquire some special calibration runs and many physics 
runs with pp collisions events at y/s = 7 TeV for an integrated luminosity of ~ 37 pb _1 . 

To illustrate the detector operation performance, we show in Fig. 2 a typical map of muon 
hits in all the five stations during a physics run. The few holes present on the maps correspond 
to dead channels. Most of them were due to hardware problems in the readout chain known since 
the beginning of the data taking 2 . They affected only 129 of the 55296 logical pads (0.2%). The 
detector failures occurred during the run, were mostly due to HV trips of single detector gaps and 
were recovered after some conditioning procedure. 

The overall effect of dead channels on the muon tracking efficiency over the full 2010 run has 
been estimated [13] by counting the fraction of muon tracks with momentum larger than 6 GeV/c 
crossing one of the dead channels in Monte Carlo minimum bias events. It resulted to be below 
1%. For the stations M2-M5, due to the large redundancy given by the four independent gaps per 
chamber, the effect of gap failures on the overall system efficiency was negligible (< 0.1%). In the 
case of Ml, gap failures are more dangerous since there are only 2 gaps per chamber. However, 
as it has been explained in section 2, this station is only used by the L0MU trigger and not by the 
high level trigger and offline reconstruction. If a detector of Ml is faulty, the L0MU trigger can be 
set to ignore the information of the trigger sector 3 containing the Ml dead channels and calculate 
Pt with M2 and M3 hits. This results in a poorer resolution of the pr measurement and a slight 
performance degradation for the affected region. This procedure was used only in one case for a 
chamber in M1R1, when both its triple GEM detectors became very inefficient. This dead zone 
would have otherwise generated an overall inefficiency of the order of 3% in the stand alone trigger 
tracking. 

The noise in the detector should be below 1 kHz per physical channel in order to have a 
negligible rate (below 1 Hz) of fake muon triggers [14]. Even lower rates are desirable to suppress 
the noise contribution to muon misidentification in the offline reconstruction. The level of detector 
noise was checked regularly with dedicated runs of random triggers in absence of beam and the 
noise rate per physical channel was computed for each detector region from the multiplicity of 
firing logical channels or logical pads. The fraction of channels having a noise rate larger than 10 
kHz is typically lower than 0.1%. Such values are well below the noise level causing a measurable 
rate of accidental triggers. 



2 Almost all of them have been cured in the LHC shutdown after the 2010 data taking. 

3 A Trigger Sector is the small detector zone which is the basic unit where L0MU operates. In the case of Ml it 
directly contains the logical pads. In most of the detector it contains horizontal and vertical strips and the fired logical 
pads are identified by crossing the x and y strips. 
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Figure 2: Illumination map of the five detector stations in a typical 2010 physics run. The log color 
scales give the average number of hits per cm 2 per trigger for all the 55296 logical pads. Faulty 
channels giving no hits and a few noisy spots can be noticed. 

5 Data samples and track reconstruction 

Due to the continous progress of the LHC, the 2010 data span a wide range of luminosity values, 
from 10 27 to 10 32 s~ ! cm~ 2 . The machine operations were accordingly evolving with time, from 
the single bunch collision per orbit of the first runs, to the 150 ns spaced bunch trains of the highest 
intensity runs with 344 bunch crossings per orbit, equivalent to a 3.9 MHz collision rate, and the 
first tests with 50 ns spaced bunches. 

In order to test the muon detector response with these changing conditions, an appropriate set 
of sample runs over the full period was chosen and analysed. 

Some special runs with events recorded in a larger time gate (125 ns, instead of the nominal 
25 ns) were also acquired for time alignment studies (time alignment events, TAE in the following). 

For the measurement of particle rates, cluster size and time resolution, only events triggered 
by some minimum bias condition, independently of the muon detector response, were used: 
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• Minimum Bias trigger (LOMB), requiring the total energy released in the HCAL to be more 
than 320 MeV; 

• "Microbias" single track trigger (microbias), requiring some hits compatible with a track in 
the VELO or first tracking stations; 

For spatial alignment studies (section 9) and detector efficiency measurements (section 10), 
when the statistics of muon tracks in the minimum bias samples was not adequate, events collected 
with standard physics triggers were used. The samples chosen and the procedures adopted to avoid 
the bias introduced by the trigger will be described there. 

Monte Carlo simulation was used to compare the observed detector performance with the 
expectation, and to verify the analysis procedures. Standard samples of LHCb Monte Carlo events 
simulating minimum bias pp interactions and production of prompt J/y decaying to pL + pi~ were 
used. Other special samples used for particular needs will be described on the next sections. The 
events were generated using PYTHIA 6.4 [15] to describe the pp collisions and GEANT4 [16] for 
the LHCb detector simulation. 

Standalone muon tracks (M-tracks) are reconstructed with an algorithm similar to the one used 
in the muon high level trigger [17]. M-tracks are reconstructed starting from the firing logical pads. 
If adjacent fired pads are found, they are clusterized, both in x and y, to obtain track hits. Clusters 
are allowed to extend across adjacent chambers even belonging to different regions. The x and y 
of the cluster barycenter and the z of the station midplane are assumed as the track hit coordinates. 
Hits aligned with the average position of the pp collision point are selected by a combinatorial 
algorithm. Track hits are fitted to a straight line and quality cuts are applied depending on the mea- 
sured quantities. M-tracks are then matched with the tracks reconstructed in the tracking detectors 
(T-tracks) when this is needed to improve the momentum resolution or reduce background. 

6 Rates 

The rate capability was one of the key request for the choice of technology and design of the muon 
detector. Detailed simulations were developed [18] to evaluate the particle rates and radiation 
doses expected for the nominal LHCb operations at an energy of y/s = 14 TeV and a luminosity 
of 2 x 10 32 s~ ! cm~ 2 . Whith respect to the radiation hardness, the detector was designed to stand 
rates larger than a factor 3 in stations M2-M5 and a factor 2 in station Ml as compared with these 
simulations, for the 10 years of planned LHCb operation. 

From the 2010 data we measured the actual rates at y/s = 7 TeV for a wide range of luminosi- 
ties and compared them with expectations [13]. From the rate seen by triggered events we could 
also extract the average number of hits per visible interaction and evaluate the contribution due to 
the current interactions and the one not due to the triggered collision. The first contribution is by 
definition independent of the luminosity. The second is expected to be dependent on the beam op- 
eration conditions. It can be due to beam background, residual detector noise or, more importantly, 
to late hits originating from collisions in previous bunches (spillover). Particle back scattering 
from the cavern or from heavy objects surrounding the beam pipe is expected to give a significant 
contribution to spillover. Also detector effects like late cross-talk or afterpulses can contribute. 

Let's consider the firing logical pads, also called hits in the following, no matter if they cluster- 
ize or not. Given an unbiased trigger T, the number of hits per unit surface and per triggered event 
rj = dNh/dSdNr was computed for every chamber. The average number of hits in each detector 
region, was calculated after removing the few chambers with channels affected by some pathology. 
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The microbias trigger provided the least biased trigger suitable for the measurement. For each 
sample, the average number /J. of interactions visible in the LHCb detector per beam crossing was 
evaluated from the fraction fo of beam crossing events not producing a trigger 



/ = P(0;ju) = e - 



(6.1) 



where P is the Poisson distribution. The pile-up factor, i.e. the average number of interactions in 
triggered events, is 

and the rate of microbias events at a luminosity Jzf is 



dN microbias -2? 

j = a x — 

dt p 



(6.3) 



where a is the cross-section evaluated to be 65 mb with 10% uncertainty, from the first luminosity 
studies. The rates for each sample are extrapolated to the nominal luminosity of 2 x 10 32 Hz/cm 2 
using the measured values of the number of firing pads per triggered event r m i cro bi as 



P ^^microbias -^nominal 

K — '"microbias X : X 



dt 



(6.4) 



The result is shown in Fig. 3. 
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Figure 3: Average hit rate, extrapolated to the nominal luminosity, in the 20 muon detector regions 
for the 2010 sample runs. The rates refer to the in-time hits (within the trigger 25 ns gate). 

The scaling of rates is verified within a few percent across five order of magnitudes in lu- 
minosity, indicating that contributions not due to the triggered collisions are small with the only 
important exception of the outer region of the last station M5 where a significant effect from back 
scattering is expected due to the limited shielding behind the detector. The contribution of back 
scattering depends on the beam conditions (jj. and bunch spacing) of each run. 

In Fig. 4 we compare the values of the average R, measured in each region for a high lumi- 
nosity run, with the predictions of the simulation. The experimental data must be compared with 
the Monte Carlo data produced in the 2010 configuration at = 7 TeV, taking in to account that 
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spillover was not simulated. The Monte Carlo reasonably reproduces the rates for the outer regions 
(except for the missing spillover contribution in M5), while underestimates them for Ml and all the 
inner regions by up to a factor 2 that is however within the safety factor assumed in the detector 
design phase. For a more complete comparison, we also report in the figure the Monte Carlo data 
produced with the fully detailed simulation including spillover, in the nominal LHCb configura- 
tion at ->/s =14 TeV. The R values are everywhere increased following the assumed increase of 
track multiplicity at higher energy and show an even higher increase in M5R4 where the spillover 
simulation reproduces a behaviour similar to the one observed in the real data. 
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Figure 4: Comparison of the rates seen in 2010 data with the MC data produced in the 2010 
configuration at y/s = 7 TeV, and the more detailed simulation, at y/s =14 TeV, used in the detector 
design phase. The rates refer to the in-time hits (within the trigger 25 ns gate). 

The spillover contribution was estimated from TAE data by using the bunch crossing identi- 
fication number (BXID) and the fine time information of the 4-bit TDC in the ODE boards. The 
LOMB events were used, since the /ibias triggers were not available in the TAE mode. A typical 
time distribution showing spillover is reported in Fig. 5. 
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Figure 5: Time distribution of the M5R2 hits for LOMB events acquired in TAE mode. The red 
vertical lines separate the consecutive 25 ns gates assigned with progressive BXID numbers. The 
structures at the gate boundaries are due to a known feature of the TDC giving an incorrect fine 
time measurement at the gate edges. 

The space distribution of any-time signals (125 ns gate) normalized to the in-time (within the 
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Contribution of late hits 




Figure 6: Space distribution of the ratio between the rates in 125 ns and 25 ns. The contribution of 
late hits from backscattering is evident in the outer regions, notably for station M5. 

standard 25 ns gate) signal distribution, is shown in Fig. 6. The plots show excess of late hits, 
clearly not related to the amount of the in-time ones, notably in the up and down edges of the more 
downstream stations. The effect is impressive in the last station M5 where the too small iron wall 
behind the detector (see Fig. 1) provides insufficient shielding from backscattered particles. It is 
worth to note that the effect of late hits is enhanced in regions where the rate of in-time hits is small 
and is suppressed where the in-time particle flux is large. This effect is visible for instance in the 
left and right edges of M5, not fully shielded from in-time particles by calorimeters and upstream 
muon stations. 

In the M5R4 region the total rate of hits in TAE events increases by a factor 10 with respect 
to the in-time rate (from 1.5 to 15 Hz/cm 2 ). The increase for the other regions can be seen in 
Fig. 7. The large effect measured for the downstream stations is essentially due to late backscattered 
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Figure 7: Dots in red: ratio of the rate in 125 ns time window to the rate in 25 ns; the out of scale 
value for M5R4 is about 10. Continuous line in blue: ratio of total cluster size to in-time cluster 
size for track hits. 



particles. The smaller increase seen for station Ml, instead, is mostly due to late cross-talk signals 
from in-time particles 4 . The overall spillover effect is somewhat large in terms of relative rate 
increase, but not worrying in terms of absolute occupancy level. It must be added that spillover 
hits, being mostly seen in the 25 ns following the bunch crossing, do not affect significantly the 
measured rates when beams are operated with 50 ns bunch crossing spacing. 
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Figure 8: Comparison of rates seen in events with beam-beam crossing (bb), single beam (eb) and 
no beams (ee). 



Finally, a check for possible sources of particles outside collision events was performed. In 
Fig. 8 the rates seen in beam-beam collision events, empty-beam and empty-empty events are 
compared. Except for the low -rate external regions of M5, rates outside collisions are less than 
10% of the collision rates. The rate of empty-beam events, triggered with the beam bunch clock, is 
due to beam gas interactions. The rate of empty-empty events, randomly triggered out of time, is 
dominated by the delayed hits from the collisions in the previous bunch. The hit rate measured far 
from LHC bunch trains is compatible with the residual detector noise. 



4 This will be clarified in section 7 were the time behaviour of cross-talk is studied. The effect of rate increase due to 
late cross-talk hits meaured for clusters associated to a track is reported in Fig. 7. 
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7 Cluster Size 



The average cluster size of muon track hits is an essential parameter of the detector response, 
since it monitors the correct operation conditions of the chambers and affects the muon trigger 
performance. In section 5 it has been described how the firing logical pads are clusterized and 
the track hits defined. The cluster size can be measured in terms of the average number of firing 
adjacent pads. Due to the finer segmentation of chambers in the bending coordinate, the average 
cluster size is sigificantly larger than 1 in x, dependent on the detector geometry of each region. A 
cross-talk between adjacent pads due to charge signal induction and to some capacitive coupling is 
expected. The cluster size values for standard operating conditions were in the past measured on 
test benches, for particles impinging perpendicularly to the chamber plane, and used to feed Monte 
Carlo simulations. 

It is worth to note that the average cluster size of muon track hits, is smaller than the one for 
all the hits in the detector, as can be seen in Fig. 9. This effect is large in regions where most of 
the radiation is due to low energy particles. If clusters associated to a track are also required to be 
isolated 5 , the cluster size is further reduced demonstrating an effect of hit coalescence. This effect 
is run dependent, being correlated with event pile-up. The isolation cut reduces the effect of pile-up 
though not fully suppressing it [13]. For this reason a run with a small /I value must be used to 
make a comparison with the cluster size measured on test benches. 
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Figure 9: Average cluster size along x in a low luminosity run for each detector region. Events are 
triggered by jiibias or LOMB triggers. Cluster selections are described in the text. 



The purely geometrical effect due to the muon trajectory inside chambers was measured by 
plotting the average x cluster size as a function of the angle between the track projection on the 
bending plane and the perpendicular to the chamber plane. A correction for this effect was per- 
formed by extrapolating to angle, as shown in Fig. 10. The resulting average cluster sizes are in 
good agreement with the design values used in the simulation as can be seen in Fig. 1 1. The cluster 
sizes never exceed 1.35, a value well within the L0 trigger requirements [14]. 

5 The isolation condition requires non-adjacent firing pads in the non-bending direction (y cluster size =1) and no 
other firing channels in the same station within 7 logical pads in x and 2 in y. 
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Figure 10: Average x cluster size for isolated track hits as a function of the track angle (in rad) for 
M2R1 (smallest logical pad region) and M5R4 (largest logical pad region). A linear fit is used to 
evaluate the cluster size for perpendicularly impinging tracks. 
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Figure 11: Average ^ cluster size at angle in experimental data and simulation (MC 2010). To 
suppress the effect of pile-up, only isolated clusters and low-luminosity data are used. 



The time behaviour of cross-talk hits was studied using TAE data. While the time distribution 
of the first pad in time in a cluster is almost fully contained in a 25 ns time window, so ensuring a 
high trigger efficiency, the other pads in the cluster can arrive significantly later as it is shown in 
Fig. 12. As a consequence, the cluster size (as the hit rate) measured in a 125 ns time window is 
larger with respect to the one mesured in the 25 ns window. The importance of the effect is region 
dependent and is quantified in Fig. 7 for clusters associated to muon tracks. The late cross-talk sig- 
nals are not normally acquired and thus do not affect the trigger performance, though contributing 
to background for future beam operations with 25 ns bunch spacing. 
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Figure 12: Time delay distribution of pads in track clusters with respect to the first pad in time: (a) 
for the small cathod pad region M1R1; (b) for the large anode pad region M5R4. 



8 Timing 



The L0MU trigger requires muon hits to be recorded in each of the five stations within the 25 ns 
LHC gate associated to a beam-beam crossing. This timing constraint is the most stringent require- 
ment for achieving the design 95% muon detection efficiency. In fact the tails in the chamber time 
response are expected to be one of the main sources of detector inefficiency. To reach the goal, the 
detectors were conceived to have a time resolution better than 4 ns at their nominal settings, while 
the 122k readout channels have to be time-aligned at the 1 ns r.m.s. level. 

The time alignment of the detector has been achieved in several steps. Test signals produced by 
a custom pulser system [10] were used for a first timing. Cosmic data collected in 2009 allowed to 
refine the channel equalization using physical signals. These two steps are described in detail in [2] 
and [19]. After the calibration with cosmic data, a satisfactory time resolution was already reached 
for all regions except the most inner ones where the statistics was a limit. The precise timing 
of the first beam particles allowed to quickly intercalibrate the channels of the highly illuminated 
inner regions. For the other regions, only a few channels exhibiting an anomalous shift in the time 
response, due to some hardware interventions, were identified and fixed. 

The detector efficiency could then be optimized by a fine tuning of the channels time offset 
with respect to the 40 MHz LHCb clock. The single channel time distribution exhibits an asymmet- 
ric shape, with a longer tail for late times. This is due to the intrinsic chamber response (dependent 
on drift time and time walk with pulse height), to the effect of delayed cross-talk hits and to the 
effect of the longer path of low momentum tracks. This implies that, in order to minimize the frac- 
tion of signals falling outside the 25 ns gate, the average time should not be centered on the middle 
of the DAQ gate, but slightly before. The optimal offset is region dependent as is the shape of the 
time spectrum. Moreover, small shifts among regions were already introduced by the optimization 
of the HV and threshold settings 

The offset optimization was then independently performed for each detector region. Special 
runs were acquired with L0MB triggered events in TAE mode, varying the global time offset in 
steps of 1 ns. For each data sample, standalone muon tracks were reconstructed requiring hits 
in all five stations. The optimal offset was chosen by maximizing the timing efficiency, defined 
as the probability that at least one of the track hits in a given station is found within the central 
25 ns gate [13]. 
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8.1 Time Performance 



The detector time performance for each detector region was estimated [13] by analysing two TAE 
event samples: 

• The first sample was acquired at the beginning of the physics data taking (before applying 
the final intercalibration for region M1R1) with events triggered by LOMB; 

• The second sample was acquired at the very end of physics data taking. At that time the 
LOMB triggers were downscaled by a factor 100 and the need of an adequate statistics re- 
quired the use of the physics L0 triggers based on calorimeters (electron, photon, hadron 
triggers). 

TAE events were fully reconstructed so that M-tracks could be required to have a good matching 
with a T-track having a momentum larger than 8 GeV/c to ensure muons to reach the M5 station. 
Ttracks were required to cross the muon detector at a safety distance from the inner and outer edges 
to avoid border effects. Clusters inside the few chambers with pathologic behaviour mentioned in 
section 4, were not considered for the timing efficiency calculation. 

The time resolution and, most importantly, the timing efficiency is evaluated from the distribu- 
tions obtained with the BXID and the 4-bit TDC mesurement for the most time centered hit in the 
track clusters. If such hit is not assigned with the correct BXID we have a timing inefficiency. As 
an example, the distribution obtained with the second TAE sample for the M5R2 region is reported 
in Fig. 13. 
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Figure 13: Time distribution of the most time centered hit measured in M5R2 for LOMB events 
acquired in TAE mode. The vertical lines delimit the "efficient" hits assigned with the BXID 
number of the trigger. 

The timing efficiencies obtained for the TAE sample acquired at the end of data taking are 
reported in Tab. 1. The overall efficiency is measured to be 98.83 ± 0.09 %, a value well beyond 
requirements. The sensitivity of the measured values to the quality cuts, suggests that the efficiency 
could be systematically underestimated, by a few per mill, in the regions most affected by combi- 
natorial background, namely region Rl and stations Ml and M5. On the other hand, the statistics 
of the available data samples does not allow to tighten the cuts, further improving the muon track 
sample purity. 

The core time resolution was evaluated by a gaussian fit to the distributions around the max- 
imum and has values between 3 and 4 ns, in line with expectations. The results obtained for the 
two TAE samples are compared in Fig. 14. Despite the different triggers used in the two samples, 
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Table 1: Timing efficiency, in per cent, for each detector region, obtained from the TAE sample 
acquired at the end of data taking. The average efficiencies on the stations are also reported. 



Station 


Rl 


R2 


R3 


R4 


Average 


Ml 


98.4 ± 0.4 


99.50 ±0.10 


99.78 ±0.10 


99.77 ± 0.26 


99.48± 0.07 


M2 


99.6 ± 0.4 


99.72 ± 0.09 


99.91 ± 0.07 


99.89 ± 0.21 


99.79± 0.05 


M3 


99.5 db 0.4 


99.72 ± 0.09 


99.91 ± 0.06 


100.00 ±0.17 


99.82± 0.04 


M4 


99.8 ± 0.4 


99.96 ± 0.06 


99.88 ± 0.07 


99.95 ±0.18 


99.91±0.03 


M5 


99.7 ± 0.5 


99.67 ± 0.10 


99.84 ± 0.07 


99.86 ± 0.20 


99.75± 0.05 



resulting in a different momentum spectrum and space distribution of the tracks, the results were 
found to be in very good agreement, except for region M1R1 where a more accurate time calibra- 
tion was used for the second sample. This demonstrates the excellent stability of the muon system 
along the 2010 run. 
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Figure 14: Core time resolution (upper plot) and timing efficiency (lower plot) measured for each 
region in the two TAE samples acquired before and after the bulk of 2010 LHCb physics data. 



8.2 Stability of the time response 

The stability of the absolute time scale in the long term is expected to be limited by two effects: 

• the LHCb clock drifts with temperature; variations were compensated manually during the 
run in order to be stable within ± 0.5 ns; 

• the variations of temperature and atmospheric pressure affect the chamber gain and variations 
of signal pulse height produce time walk effects. The largest effect is expected from pressure 
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variations, and is estimated to be equivalent to a ~ ± 20 V change in HV [20], corresponding 
to ~ ± 0.4 ns. 



The average time of track hits, measured as a function of time in the sample runs acquired across 
the 2010 run, is shown in Fig. 15 for various detector regions. The behaviour is consistent with the 
mentioned effects. Variations are at the level of ± 1 ns and are clearly correlated among regions. 
Residual uncorrelated variations are compatible with zero. There is no evidence for a dependence 
of the time drift on the detector illumination, that could be a hint for an ageing effect on the chamber 
gain. 
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Figure 15: Variation of muon chambers time response during the 2010 run. The average time of 
the most time-centered track hit, in different stations and regions, is plotted as a function of the 
data-taking time. The values for optimal efficiency depend on the detector region and have been 
fixed at the start of data taking. Variations along the 7 months of operation do not exceed the ± 1 ns 
range. 



9 Spatial alignment 

The spatial alignment of the muon detector must guarantee the design performance of trigger and 
offline muon identification. As explained in section 2, the L0MU trigger requires hits in all the 
5 stations aligned on a muon track-segment having a pj above a given threshold. Offline muon 
identification is more flexible in the reconstruction of the muon track-segment but requires the 
matching with a track reconstructed in the tracking system through the whole spectrometer. 

The alignment accuracy needed in the system is driven by the trigger requirements in the fine 
granularity inner regions of the stations Ml, M2 and M3. In particular in these three stations the 
trigger FOI is defined by 1 pad only in the non-bending vertical coordinate y and a misalignment can 
directly contribute to trigger inefficiency. In the bending coordinate x, the FOIs are composed of 
several pads and the main effect of a misalignment can be a bias in the pj calculation if the trigger 
algorithm does not use the true x positions. The alignment of stations M4 and M5 is less important 
because in their case the y FOI is as large as 3 pads and their hits are not used to calculate pj. The 
detector mechanics was designed to reach a precision of the order of 1 mm in x and y directions. 
The alignment requirements along z are much less demanding due to the forward geometry of the 
experiment. 
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During the installation, the supporting walls were kept in the open position and the muon 
chambers mounted with an accuracy of ~ 1 mm along x and y coordinates, centered on their nominal 
positions, relative to reference targets placed on top of each half station. The measured rotations 
were zero within the precision of 1 mrad. After chamber installation, the half stations were closed 
around the beam pipe leaving a small safety distance between the A and C side; the two half 
stations being ideally positioned left-right symmetric and projective to the interaction point. The 
closed stations were then surveyed with respect to the LHCb cavern using four reference targets 
on each side, and the values stored in the geometry database used by the offline reconstruction 
program to define the absolute hit coordinates. The values measured by the survey for the 2010 
collision run are reported in Fig. 16. The MIC half station shows a non negligible though small 
misalignment from the ideal position. 
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Figure 16: Alignments of the ten muon half stations for the 2010 run. The average value x of the 
inner edges (left) and the median y (right), are shown as a function of the station z position. The 
empty dots represent the survey measurements, the full dots are the positions found by the software 
global alignment. The error bars correspond to the statistical and systematic uncertainties, summed 
in quadrature. The dashed lines in the left plot represent the ideal alignment. 



9.1 Space alignment measured with muon tracks 

A study of the muon system alignment was performed using muon tracks with the first purpose of 
checking the mechanical positioning and the survey measurements. In the same time a tool was 
prepared for alignment monitoring after each intervention requiring the opening of the stations, and 
eventually for correcting the geometry data base. 

In this analysis, standalone muon track segments (M-tracks) are defined by at least four clus- 
ters in four different stations that are compatible with a straight line. Candidates with too large 
clusters (number of pads>6) are eliminated to avoid mis-reconstruction problems. Moreover sta- 
tions with more than 300 hits are excluded. The M-tracks are required to match a good quality 
track reconstructed in all the tracking detectors (long T- track) having a momentum p >6 GeV/c. 
The matching condition requires a good % 2 between the parameters of M-track and T-rrack extrap- 
olated to the M2 position. Matching segments are then merged together in a unique track that is 
required to have a good % 2 . 

The alignment is performed following the standard procedure used in LHCb [21]. It makes 
use of the Kalmanfit iterative method, that performs a minimization of the total % 2 of an ensemble 
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of tracks while adjusting the detector positions. Only the positions of the muon half stations were 
allowed to vary since the tracking detectors were previously independently aligned using the same 
procedure. The iterative process starts assuming the muon detector in the positions measured by the 
survey and stops when the convergence is reached (total % 2 doesn't improve significantly) usually 
after 4-5 iterations. In Tab. 2 are reported the misalignments relative to the survey, measured on 
a sample of ~ 7000 tracks, with an analysis performed fitting only the translational degrees of 
freedom along x and y of each half station. The absolute positions in the LHCb reference system 
are shown in Fig. 16 together with the survey measurements. 

Table 2: Misalignments of muon half stations M1-M5, relative to the survey measurements, calcu- 
lated in the LHCb reference system with the Kalman fit iterative method. The quoted uncertainties 
are the fit errors (first) and the systematic uncertainties (second) on the relative positions determined 
repeating the analysis with different track selections. 



Ml 


C-side 

Ax (mm) Ay (mm) 
0.92±0.14±0.15 0.11±0.28±0.24 


A-side 

Ax (mm) Ay (mm) 
1.04±0.14±0.24 0.68±0.27±0.20 


M2 


-1.56±0.05±0.04 -0.89±0.12±0.13 


-0.42±0.05±0.07 0.76±0.12±0.17 


M3 


-2.41±0.08±0.05 -1.53±0.14±0.14 


-1.19±0.08±0.06 0.45±0.14±0.08 


M4 


-0.33±0.15±0.11 -1.59±0.17±0.09 


-1.44±0.14±0.20 0.28±0.17±0.15 


M5 


-2.14±0.18±0.12 0.18±0.20±0.11 


-0.07±0.18±0.20 0.16±0.20±0.15 



The systematic errors have been estimated repeating the analysis with samples of tracks hav- 
ing different momenta, comparing subsamples of tracks hitting different detector regions, and using 
slightly different alignments for the tracking system. While the absolute positions of the stations 
show variations of the order of 1 mm both in x and y, their relative positions are more stable. The 
results obtained with the Kalman fit iterative method show significant even though small misalign- 
ments relative to the survey measurements (in particular for the x position of Ml, M2, M3 side C) 
and suggest to apply a correction to the geometry data base. Studies on the additional degrees of 
freedom like rotations around y or x directions give results compatible with the survey. Due to the 
forward geometry of the experimental apparatus, shifts in z direction are as difficult to measure as 
unimportant, therefore the z values measured by the survey were assumed. 

10 Detector efficiency 

The overall performance of the muon detector is quantified by the detection efficiency of muon 
tracks when the system is operated in the standard data taking conditions. The inefficiency intro- 
duced by dead channels, or other hardware failures occurred in the 2010 detector operation, has 
been quantified in section 3. Here the intrinsic efficiency of the system is evaluated after applying 
strict fiducial volume cuts and eliminating the few small zones where known problems are present. 
The used procedure is described in detail in reference [13] and is here summarized: 

• Different data sets are used for Ml and M2-M5 stations to select appropriate samples of trig- 
ger unbiased standalone muon tracks (M-tracks). To reach the needed purity of the sample, 
M-tracks are matched with a good quality, high momentum T-track and required to fulfill 
tight selection criteria (section 10.1). 
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• The presence of background hits affects the results and requires different procedures to cor- 
rectly evaluate the true efficiency for M2-M5 stations and for M 1 station where the occu- 
pancy is much higher (section 10.2). 

• The efficiency for each station is estimated by searching clusters around the prediction de- 
fined by the M-track reconstructed using only the other 4 stations. The search of clusters 
around the prediction is repeated increasing the opening window from 1 up to 8 standard 
deviations, both in x and y. The value of a being determined, region by region, by a gaussian 
fit to the central part of the distribution of the distance between the position predicted by the 
M-track and all clusters in that region. For Ml the prediction is defined by the T- track asso- 
ciated to the M-track in order to improve its quality. The values of a for the twenty regions 
are reported in Tab. 3. 

Table 3: Resolution along x and y of the distance between the muon track and the muon cluster 
in each region of the muon detector. The muon track is reconstructed skipping the station whose 
resolution must be evaluated. 







Ml 


M2 


M3 


M4 


M5 


Rl 


o x x a v (mm 2 ) 


4x 10 


15x30 


10 x 12 


15 x 16 


33x40 


R2 


o x x o y (mm 2 ) 


8x 18 


25x50 


15 x24 


27x32 


50x60 


R3 


a x x o v (mm 2 ) 


16x40 


35 x70 


25x48 


48 x 64 


100 x 110 


R4 


a x x Oy (mm 2 ) 


32x80 


60 x 100 


40x96 


97 x 128 


150 x 180 



If the full procedure is correct, the efficiency measured as a function of the search opening 
window will show a saturation behaviour permitting a reliable estimate of the detector efficiency. 
Particular care is required in the background subtraction that must be a stable and reliable procedure 
even when a search area as large as 16 <7 A x 16 o y is considered. 

10.1 Muon samples and track selection 

Different data samples are used for Ml and for M2-M5 stations: 

• For the efficiency of M2-M5 stations, data acquired in two fills, corresponding to an inte- 
grated luminosity of 1.2 nb _1 , were used. The first level hardware trigger (L0) required a 
high pj hadron or lepton detected in the calorimeter or in the muon system; the first step 
of the software high level trigger (HLT1) required the logical "or" of several independent 
algorithms. To remove the bias introduced by the trigger in the efficiency calculation, events 
where both L0 and HLT1 were fired irrespectively of the muon system information were 
selected. With this data sample, the majority of the muons reaching the muon stations and 
used for the analysis, originate from decays in flight of 7r's or K's. 

• For the Ml efficiency measurement, kaons decaying at the end of the tracking system can 
generate a good T-track giving a poor quality Ml prediction, not adequate to the large 
occupancy of the station. To have a sample of true muons, events with a reconstructed 
J/ljf—^ M + M were used. This sample corresponds to almost all data acquired in 2010 
(~37pb- ! ). To remove the L0 bias on the efficiency evaluation, in each li + pL pair from 
a. J/y, the muon which fired the L0 trigger was not considered in the analysis. Notice that 
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the use of the J/xjf sample for the analysis of M2-M5 station is not possible since the muon 
information is used to identify both muons in the HLT and in the reconstruction. 

A tight selection is required to reach the purity of the M-track sample needed for a precise efficiency 
measurement and different conditions are required for the different stations. 

For every station, when its efficiency is being evaluated, the M-track is validated by requiring 
the matching with a T-track in the other four stations used for the fit. The matching requires that 
the distance from the T-track extrapolated to the station and the clusters associated to the M-track 
is within one standard deviation in both x and y projections. Such standard deviations are estimated 
by a gaussian fit to the central peak of the distance distributions. 

A momentum cut of 12 GeV/c (15 GeV/c) is applied to the T-track when M2 (M3, M4, M5) 
station is analysed. If more than one T-track matches the M-track candidate and at least one of 
them has a momentum below the cut, the candidate is rejected. When evaluating the efficiency of 
M2 and M3, where the occupancy is relatively high and the fired hits are identified by crossing 
vertical and horizontal strips, further cuts on the local hit multiplicity are applied to avoid ghost 
combinations. Moreover, when analyzing the efficiency of M3 where the prediction resolution is 
poor due to the lower granularity of M4, a cluster size of 1 is required on M2 station. For the 
analysis of Ml station, the M-track sample selection starts from the T-track associated to a muon 
candidate from the J/y and a cut on momentum of 12 GeV/c is applied. 

10.2 Background subtraction 

The presence of background clusters affects the search results. The necessary background subtrac- 
tion requires different procedures for M2-M5 and for M 1 station where the occupancy is higher. 

Assuming a Poissonian nature of the background, the average number of clusters due to the 
background and the efficiency in a given search window can be extracted by a fit to the distribution 
of the number of clusters found. If this method works, the background estimate in the neighbouring 
of the muon track takes automatically into account any possible correlation between muons and 
background, as in the case of delta rays or punch through in the calorimeter. 

For stations M2-M5, a fully satisfactory result is obtained assuming a background with two 
Poissonian components; the probability of finding n clusters in the search window being 

^^l^+^'i^i+^H-^+u-)-^] 00.1) 

where e is the efficiency, B 1 and B2 are the two Poissonian components of the background and r 
their ratio. 

As an example, in Fig. 17 are reported the results of the fit in the case of an opening window 
of 8a r and 8cy for the station M3. It is worth to add that a fit with only one Poissonian component 
for the background does not give an equally good representation of the high multiplicity bins, but 
the fitted values of the efficiency do not show significant differences. 

In the case of the more crowded station Ml where the J/y muon sample is used, the cluster 
multiplicities are not very well fitted assuming a background with two Poissonian components. 
Thus another fully independent method is also used to evaluate the background and the efficiency. 
It exploits the <p rotation invariance of the primary interactions and assumes that the background 
correlated with the muon track is negligible. The background is estimated by counting the average 
number of clusters in the same search window but in the opposite quadrant, with respect to the 
track prediction. The soundness of the method is confirmed by comparing the number of clusters 
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Figure 17: Multiplicity of clusters found in the 8(7 search window for the M3 regions. The thick line shows 
the results of the fit with the two components background described in equation 10.1. 



found in the opposite quadrant with the number of clusters in the track prediction quadrant having 
subtracted one cluster attributed to the muon track, as shown in Fig. 18. 
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Figure 18: The cluster multiplicity in the 8c search window diminished by 1 (histogram) and the cluster 
multiplicity in the corresponding window in the opposite quadrant (dots), for J/y/ muons in the four regions 
of Ml station. 



The true efficiency e t is then estimated by the formula: 

e a = e t + (l-e t )-P bg (10.2) 
where e a is the apparent efficiency calculated as NNcws>o/Np re d s , being Nnclus>q the number of 



tracks where at least one cluster has been found and Np re d s the total number of tracks predicted to 
fall in the search window; Pj, g = Nncws>o,oqI Np re d s is the probability to find at least one cluster 
in the search window in the opposite quadrant. 

The efficiency values obtained with this method do not show any significant difference with 
the values obtained with the previous one. 

10.3 Check of the procedure with simulation 

To evaluate if the whole procedure of track selection, background subtraction and efficiency deter- 
mination, is robust and unbiased, a test has been made using events simulated with a Monte Carlo 
where all background components and detector effects are described. 

The procedure was applied to the Minimum Bias (for M2-M5) and J '/y (for Ml) Monte 
Carlo events and the efficiencies compared with the ones obtained applying the same procedure to 
an ideal sample of muon tracks the so-called Particle Gun (PG) muons. These muons are generated 
as starting from the interaction point with a predefined momentum and angular distribution. Since 
those events are practically background free and do not suffer from fakes in track reconstruction 
and selection, they allow to cleanly extract the detector efficiency with the present method. Inde- 
pendently of the absolute values found for the two samples, a satisfactory agrement would then be 
an overall check of the correctness of the procedure. The agreement was found to be good for the 
stations M1-M4, while for R3 and R4 regions of the M5 station the efficiency extracted from the 
Minimum Bias is lower than the corresponding PG efficiency as shown in Tab. 4. This difference 
is explained with a lack of purity of the Minimum Bias track sample: part of the selected tracks 
are muons produced in the showers inside the calorimeter that are aligned and matched with the 
hadron track initiating the shower. These muons have sometimes low momentum and are absorbed 
between M4 and M5. This effect is region dependent 6 and can cause an artificial inefficiency in 
M5. Removing the muons with momentum below 6 GeV/c, the efficiencies reach values in agree- 
ment with the PG sample as can be seen in Tab. 4. To take into account such effect, the real data 
results will be corrected for the ratio between the efficiencies with and without the 6 GeV/c cut, as 
extracted by the Minimum Bias Monte Carlo sample. 

Table 4: M5 Monte Carlo efficiency (%) for Minimum Bias (MB) events with and without 6 GeV/c 
momentum cut compared to Particle Gun (PG) muons. 





Rl 


R2 


R3 


R4 


PG 
MB 

MB (p > 6 GeV/c) 


99.5 ± 0.3 
99.58 ± 0.07 
99.77 ±0.07 


99.73 ± 0.09 
99.46 ± 0.04 
99.63 ±0.04 


99.69 ± 0.06 
99.14 ±0.07 
99.51 ±0.08 


99.88 ± 0.03 
99.0 ± 0.2 
99.8±0.1 



10.4 Measured efficiencies 

The behaviour of the efficiency as a function of the number of d's of the opening window has been 
analysed for each region of the muon system. In all cases a correct saturation is observed at 3-4a 
demonstrating the reliability of the method. However the final value taken for the efficiency is the 
one at 8 a to allow for the presence of non gaussian tails in the prediction point. 

6 Tracks crossing outer regions have, on average, lower momentum than those crossing inner regions. 
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In Fig. 19 is shown the behaviour of the efficiency as a function of the number of a for the 
four regions of M2 station. In Fig. 20 the same is shown for the station Ml, with the efficiencies 
obtained with equation 10.2. As anticipated, the efficiencies obtained for Ml with the other method, 
assuming a background with two Poissonian components, are in excellent agreement. The average 
of the two estimates was assumed for the efficiency. 
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Figure 19: The measured efficiency of M2 station as a function of the number of a of the search window. 
The four regions are shown: Rl (red circles), R2 (green squares), R3 (black triangles), R4 (pink crosses). 
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Figure 20: The measured efficiency of Ml station as a function of the number of a of the search window. 
The four regions are shown: Rl (red circles), R2 (green squares), R3 (black triangles), R4 (pink crosses). 



In Tab. 5 the measured efficiencies with the statistical and systematic errors are reported. For 
M2-M5 the systematic errors due to background modeling have been estimated by changing the fit 
function from two poissonians to a single one. The efficiency varies less than 0.01% in all regions 
of M2-M5. For Ml where two independent methods of estimating the background have been used, 
a systematic uncertainty of half the difference between the two results has been assumed. Since 
the choice of evaluating the final efficiency at 8a has a certain degree of arbitrariness, a systematic 
error of half the difference between the efficiency value calculated at 4a and 8 a has been assumed. 
The uncertainty, due to MC limited statistics, of the correction applied on the M5 efficiency to take 
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into account the absorption of muons between M4 and M5, is included in the systematics. The 
different systematic sources have been added in quadrature. 



Table 5: Efficiency values (%) of all the twenty regions with their statistical and systematic errors. 





Rl 


R2 


R3 


R4 


Ml 


98.66 ± 0.08 ± 0.08 


99.37 ± 0.04 ±0.08 


99.70 ± 0.03 ±0.06 


99.85 ± 0.04±0.02 


M2 


99.68 ± 0.05 ±0.04 


99.78 ± 0.03 ±0.05 


99.79 ± 0.04 ±0.02 


99.8 ±0.1±0.01 


M3 


99.35 ± 0.06 ±0.05 


99.79 ± 0.02 ±0.04 


99.85 ± 0.02 ±0.01 


99.86 ± 0.06 ±0.01 


M4 


99.62 ± 0.03 ±0.07 


99.89 ± 0.01 ±0.03 


99.65 ± 0.04 ±0.03 


99.72 ± 0.09 ±0.01 


M5 


99.64 ±0.07 ±0.10 


99.82 ± 0.03 ±0.04 


99.90 ±0.04 ±0.12 


100.0-0.2-0.1 



The overall detector efficiency is compared in Fig. 21 with the results of section 8.1 where the 
contributions of timing tails to the inefficiency are evaluated. As already mentioned, the limited 
statistics of TAE data prevented to apply the same track selection for the two studies. The lower 
purity of the muon track sample for the timing analysis leads to a slight underestimate of the 
timing efficiency. The data were also acquired at different times, and this could explain the poorer 
performance of region M1R1 for the TAE data acquired at the end of the runs with a lower high 
voltage (see section 4). Nevertheless, the region dependence of the two measurements is very 
similar, and for regions less affected by combinatorial background the timing efficiency is found 
to be compatible or only slightly better than the overall efficiency, indicating that the bulk of the 
detector inefficiency comes from signals falling outside the 25 ns LHC gate. 
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Figure 21: The results for total efficiency are compared with the estimates of timing efficiency from 
section 8.1. Errors are statistical only. For the regions more affected by combinatorial background, 
a systematic underestimation of the timing efficiency by a few per mill is present. 
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11 Conclusions 

The muon detector was successfully operated since the first year of LHC physics. Its performance 
has been studied across the five orders of magnitude of luminosity experienced during 2010 and 
compared with expectation. The whole system demonstrated an excellent reliability and stability. 
Detector requirements in terms of cluster size, time resolution and efficiency are fullfilled. 
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During 2010 a small number of dead channels accounted for an overall inefficiency in the 
muon detection < 1%. Most of these channels were cured in the LHC shutdown after the 2010 
data taking. Thanks to a good monitoring system and maintenance work, temporary failures oc- 
curring during the run, mainly HV trips of single chamber gaps, accounted for per mil effects in 
the efficiency. A careful setting of the chamber working point and a precise timing intercalibration 
allowed to reach a muon detection efficiency, mainly determined by the chamber time resolution, 
well above the design requirement of 99% in all the 5 muon stations. 
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